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The measurements of the centrality dependence of dN/dr) and transverse momentum spectra 
from mid- to forward rapidity in d+Au collisions at ^fs~^= 200 GeV are reported. They provide 
a sensitive tool for understanding the dynamics of multi-particle production in the high parton- 
density regime. In particular, we observe strong suppression of the nuclear modification factor Rcp 
at forward rapidities (d-side, n — 3.1) and enhancement at backward rapidity (n — —3.1). An 
empirical scaling is obtained for multiplicity and Rcp when a shift of the center-of-mass in the 
asymmetric d+Au collisions with respect to the nucleon-nucleon system is applied. 

PACS numbers: 25.75.-q, 25.75.Dw, 13.85.-t 



I. INTRODUCTION compared to Au+Au collisions. Measurements of the nu- 

The d+Au collisions at the Rclativistic Heavy Ion Col- 
lider (RHIC) provide an important control environment 
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clear modification factor and back-to-back correlations 
at mid-rapidity in d+Au collisions suggest that the sup- 
pression of particles with high transverse momentum and 
the disappearance of back-to-back correlations, seen in 
Au+Au collisions are due to final-state interactions 
with the hot, dense medium produced in such collisions, 
rather than initial-state effects on the Au nucleus @ ■ The 
observed enhancement of the nuclear modification factor 
in the region of transverse momentum ivt) of 2 GeV/c 
at mid-rapidity in d+Au collisions 0, 0, la, E| j referred 
to as the "Cronin effect" Q, can be described within a 
pQCD framework incorporating initial multiple parton 
scattering and nuclear shadowing d, Q . Saturation ef- 
fects (mostly described as the formation of Color Glass 
Condensate (CGC) are expected to be more pro- 

nounced at large rapidity (y) or pseudorapidity (77) close 
to the deuteron beam, where the small- 2; components of 
the Au nucleus wave function can be probed. Recent re- 
sults reported by the BRAHMS collaboration [ll[ , where 
a suppression of the nuclear modification factor at for- 
ward rapidities is visible, are in qualitative agreement 
with predictions within the framework of gluon satura- 
tion in the CGC [l2|. These results indicate a possible 
dramatic evolution of gluon saturation from mid- to for- 
ward rapidity at RHIC. However, it should be noted that 
these results can be reasonably described by pQCD mod- 
els [1, H| and in the framework of final-state parton re- 
combination (l3j . On the other hand, as the rapidity of 
the probe decreases (and at the same time, x increases), 
the multiple-scattering contribution to the Cronin effect 
should decrease. The first results from the PHENIX col- 
laboration [Tij show an opposite behaviour which cannot 
be explained by current model calculations [l5l |. 

In this paper the pseudorapidity and centrality depen- 
dence of the nuclear modification factor Rcp will be dis- 
cussed in connection with the asymmetry in particle pro- 
duction in d+Au collisions at ^/snn =200 GeV. 



II. EXPERIMENTAL SETUP 

The STAR experiment [l6[ at RHIC measures charged 
hadrons over a wide range of pseudorapidity and trans- 
verse momentum. The main detector is a large Time Pro- 
jection Chamber (TPC) [T3| which allows particle identi- 
fication via dE/dx within the range of \rj\ < 1. Charged 
particle detection in the forward directions is achieved 
with the two azimuthally symmetric Forward TPC (FT- 
PCs [lH) which extend the pseudorapidity coverage of 
STAR to the region 2.5 < \tj\ < 4. The FTPCs which 
utilize a radial drift field perpendicular to the magnetic 
field, achieve a two-track resolution of 2 - 2.5 mm (an 
order of magnitude better than a TPC using a constant 
drift field). This allows track reconstruction in the larger 
rapidity region where track densities are high. In d+Au 
collisions at a ccntcr-of-mass energy of y/sNN=200 GeV 
a track finding efficiency in the FTPCs of about 90% was 
reached, independent of centrality, 77 and pt in the phase 



space region of 3 < 77 < 3.5 and 0.1 < p t < 3 GeV/c used 
for this analysis. The momentum resolution in the FT- 
PCs is also independent of centrality, but shows a strong 
dependence on 77 andpt- For T) « 3.1, the relative momen- 
tum resolution degrades approximately linearly from 10% 
to 25% - 30% in the region of 0.1 < p t < 3 GeV/c and 
the pseudorapidity resolution is better than 0.02 units 
in 77. Background and secondary decay products cor- 
rections were estimated using HIJING simulations [lj|. 
The main systematic error quoted in this analysis (if not 
otherwise mentioned) is caused by the momentum resolu- 
tion of the FTPCs. This affected mainly the vertex DCA 
(distance of closest approach) requirement used to select 
primary charged hadrons in the FTPCs. An estimate of 
the main systematic error was done by varying the DCA 
by ± 0.5 cm. A detailed description of the various cal- 
ibration steps, further corrections and data quality can 
be found in [201 ] . 



III. MEASUREMENTS 

A. Mean transverse momentum 

To quantify the influence of momentum resolution on 
the transverse momentum spectra simulated charged pi- 
ons were embedded in real d+Au events, adding up to 5% 
of the total event multiplicity in the FTPCs. Initially flat 
input distributions were weighted according to the mea- 
sured transverse momentum distributions as function of 
centrality and transverse momentum (20| . Background 
and secondary decay products corrections were estimated 
using HIJING simulations 0, [2(| . The transverse mo- 
mentum distributions corrected for momentum resolu- 
tion, background and secondary decay products mea- 
sured in the East-FTPC (Au-side; 77 « -3.1) and West- 
FTPC (d-side; 77 w 3.1) arc shown in Fig. [TJi for dif- 
ferent centrality classes. Fig. [Tp shows the mean trans- 
verse momentum (p t ) obtained from a power-law fit to 
the measured p t spectra in the region of 0.1 < pt < 1 
GeV/c as a function of the number N part of participant 
nucleons. In contrast to the "naive" picture in which the 
partons from the d-nucleus experience multiple collisions 
while traversing the Au-nucleus, and therefore ac quir e 
an enhancement of (pt) on the d-side of the collision [Tj|, 
a slight increase of (pt) with centrality is visible on the 
Au-side, whereas on the d-side virtually no centrality de- 
pendence of (p t ) is present. 



B. Rcp in the forward directions 

Another variable of interest is the ratio of central to 
peripheral inclusive d+Au spectra 

ff , {d 2 N/dp t drj/{N bi n))\central 

n CP (Pt) - {d 2 N/dptdrj/{Nun) )i per . ph > y > 
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FIG. 1: (Color online) a) Charged hadron transverse momentum distributions at |r/| ~ 3.1 for 0-20%, 20-40% and 40-100% 
central d+Au events. Closed triangles refer to measurements on the Au-side (rj m —3.1) and open triangles to the d-side (rj w 
3.1). b) (pt) extracted from fitting the transverse momentum spectra with a power law function at \r]\ ~ 3.1 as a function 
of Npart- The closed triangles are measurements on the Au-side and the open triangles on the d-side of a d+Au collision. 
Statistical (lines) and systematical (bars) are shown separately (for details see text). 
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FIG. 2: (Color online) Rcp at rj m —3.1 (Au-side; triangles) and at i) ~ 3.1 (d-side; circles) for 20% and 20-40% central d+Au 
events. The grey band represents the estimate of the systematic error [2Cl |. Midra pid ity charged hadron Rcp (0-20%) (|?y| < 
0.5) is shown as stars [U. Measurements from BRAHMS at rj w 3.2 on the d-side [ill ] and from PHENIX in the region of —2 
< 77 < —1.4 on the Au-side [3| are also overlayed. The solid lines represent HIJING simulations [lj| without shadowing, the 
dashed lines with shadowing. 



where d 2 N / dptdrj is the differential yield per event and 
(Nbin) the mean number of binary collisions for the corre- 
sponding centrality class, calculated using a Monte Carlo 
Glauber model Q (see Table HJ). For the Rcp mea- 
surements the transverse momentum range could be ex- 
panded to 3 GeV/c (instead of 1 GeV/c for the trans- 
verse momentum measurements in section IIII A[) due to 
the centrality independence of the momentum resolu- 
tion. Therefore, the effect of momentum resolution on 
the transverse momentum spectra cancels out in the Rc p 
ratio. 



Centrality class 


■A^ part 


Nbin 


0-20% 


15.69 ± 1.19 


15.07 ± 1.29 


20-40% 


11.17 ± 1.11 


10.61 ± 0.80 


40-100% 


5.14 ± 0.44 


4.21 ± 0.49 


0-100% 


8.31 ± 0.37 


7.51 ± 0.39 



TABLE I: N par t and Nbi n for various centrality classes. 



Comparing forward to backward rapidities, it can be 
seen in Fig. [5] that Rcp is increasing with p t for p t < 3 
GeV/c on the Au-side. Also Rcp is larger on the Au-side, 
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FIG. 3: (Color online) Charged hadron pseudorapidity distribution per event in the TPC and FTPC acceptance for 0-20%, 
20-40%, 40-100% central and minimum bias d+Au events (triangles). The error bars include both statistical and systematic 
error. Also overlayed are measurements from BRAHMS (crosses) [22| and PHOBOS (circles) [23 ;. 24]. In addition predictions 
from HIJING AMPT [25| and the saturation model [l^] are also plotted compared with STAR measurements. 



which indicates that the Cronin effect is more pronounced 
on the Au-side of a d+Au collision. Strong centrality de- 
pendence of Rc p on the d-side of the collisions is another 
interesting feature seen in Fig. [21 This observation was 
also reported by the BRAHMS collaboration . How- 
ever, no significant centrality dependence is observed on 
the Au-side of a d+Au collision. The Rcp measurements 
from the FTPCs are in good overall agreement with mea- 
surements from BRAHMS on the d-side [ll| for p t > 1 
GeV/c and in agreement with PHENIX on the Au-side 
[lH of a d+Au collision (see Fig. The discrepancy 
between the BRAHMS and STAR Rcp measurements 
on the d-side at t] f» 3 for p t < 1 GeV/c can not be 
completely resolved at this time. The discrepancy can be 
partially attributed to the different centrality classes used 
by the two experiments to calculate the Rcp- BRAHMS 
uses 60-80% as the most peripheral bin whereas STAR 
uses the 40-100% centrality class. Also, the BRAHMS 
centrality selection is biased towards peripheral collisions 
in forward rapidities as discussed in section [III CI Fur- 
thermore, different low-pt cut-offs may affect the low pt 
measurements, where the difference between BRAHMS 
and STAR is most prominent. 

The suppression of Rcp (and RdAu [HJ) at higher ra- 
pidities on the d-side is in qualitative agreement with pre- 
dictions of the saturation model [l^ . Models based on 
pQCD which incorporate initial-state parton scattering 
and energy loss can also describe the behaviour of Rcp 
at higher rapidities [8, 9]. Furthermore, in the frame- 
work of parton recombination in the final state, Rcp at 
forward rapidities can be described as well fl3j . 

C. Charged particle density asymmetry 

In Fig. [3^, the pseudorapidity distribution dN c h/drj of 
charged hadrons per event in the TPC and FTPC accep- 
tance is shown for minimum bias and for the 0-20%, 20- 



40%, and 40-100% most central events. For comparison, 
measurements of the pseudorapidity distribution from 
BRAHMS [12 and PHOBOS [^l, S5[ are also plotted. 
The measured dN c h/dr] distributions for minimum bias 
d+Au events are in good agreement for all three experi- 
ments. However, with increasing centrality, a significant 
difference in the particle density at negative pseudora- 
pidity values 7] < -3 between STAR and PHOBOS is 
visible. On the other hand, the measurements in the 
mid-pseudorapidity region are in good agreement. When 
comparing central events, for the r/ < —3 region the 
BRAHMS dN ch /dri distribution is lower than the STAR 
measurements; at mid-rapidity it is higher. A possible 
explanation could be the different methods used for cen- 
trality selection. Centrality selection for the STAR-TPC 
was done via the N c h multiplicity in the FTPC and vice 
versa, to avoid autocorrelations caused by fluctuations in 
the measured multiplicity. Simulation studies show that 
with a pseudorapidity gap of 2 units between the detec- 
tors this method is insensitive to autocorrelations [2fj| . 
Use of the FTPC N c h multiplicity on the Au-side instead 
leads to a visibly higher particle density in the dN c h/dr] 
distribution for r\ < —3. causing a significant bias in the 
centrality definition [2fj|. This observation explains the 
higher particle density measured by PHOBOS for the 
Au-side of a d+Au collision, because their centrality was 
determined via the multiplicity in the pseudorapidity re- 
gion of —4 < rj < —3.5. For BRAHMS the enhancement 
in the particle density at midrapidity could be due to 
the fact that the multiplicity in the central region \r)\ < 
2.2 was used to define centrality. However, within the 
systematic errors the results of all three experiments are 
consistent with each other. 

In addition, the measured pseudorapidity distributions 
were compared with model predictions. Calculations 
based on gluon saturation in the Color Glass Condensate 
[3 as well as results of HIJING [H and a Multi-Phase 
Transport Model (AMPT) HH are shown in Fig. \5]p. All 
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model calculations are in good overall agreement with the 
measured dN^/drj distributions for different centrality 
classes. In particular, the models are able to reproduce 
the increasing asymmetry of charged particle densities 
with increasing centrality. 



IV. COLLISION ASYMMETRY AND YIELD 
SUPPRESSION 

In the following section the centrality dependence of 
Rcp at high rapidities |?7| « 3.1 will be discussed in con- 
nection with the observed increasing asymmetry of the 
produced particle density in d+Au collisions at trans- 
verse momenta p t < 3 GcV/c. In Fig. [2 the Rcp from 
HIJING simulations for different centralities with and 
without shadowing is shown. It is evident from Fig. [5] 
that HIJING reproduces the overall behaviour of Rcp 
at \rj\ « 3.1. In addition, it can be concluded that the 
influence of shadowing for pt > 1 GeV/c only affects the 
measurements on the Au-sidc. Since the pt spectrum 
on the d-side is more or less independent of centrality 
- except for an overall scale - (see Fig. QJi and b) and 
comparable with p+p collisions, the suppression of Rcp 
on the d-side (see Fig. [5]) could be due to the asymmetry 
in particle production in d+Au collisions with respect to 
the symmetric p+p collisions (see Fig. HJ. To take the 
asymmetry in d+Au collisions into account, a new vari- 
able t]cm is introduced, which is defined as the weighted 
mean of the dN c h/drj distribution for each centrality class 
(i?CM = VCM{N part )). r/cM was extracted from the pub- 
lished PHOBOS results [24| and should represent the 
shift of the center-of-mass in the asymmetric d+Au colli- 
sions with respect to the nucleon-nucleon center-of-mass 
system (see Tabic [TTJ) . Even though the PHOBOS data 
are biased towards higher multiplicity on the Au-side as 
discussed in section MI CI they were used to determine 
Vcm to maintain a model independent approach (also 
they are the only measurements available covering the 
full T) ran ge) . The dN c h/dr] distribution for inelastic p+p 
collisions 



A 3.! 
CM 



22j in this new reference system - obtained by 
shifting with rjcu - are shown in Fig. [4j One observes 
that the d+Au dN c h/drj distributions normalized with 
(Npart/2) (see Table Q} at high rapidities are consistent 
with the shifted p+p values. Also, the centrality depen- 
dence can be qualitatively explained. Therefore, r]cM 
seems to be an appropriate variable to describe the asym- 
metry in particle production in d+Au collisions assuming 
that this asymmetry is caused by the nuclear stopping of 
the deutcron while traversing through the gold nucleus. 
Similar approaches to describe the pseudorapidity distri- 
butions in d+Au can be found in [2(| H3, Ha] . 

In that representation, the suppression of the parti- 
cle density on the d-side and the enhancement on the 
Au-side in asymmetric d+Au collisons with respect to 
the symmetric p+p collisions (see Fig. @| can be ex- 
pressed by defining SdAuiViVCAi) as the ratio of the 
p+p reference dN^/drj distribution shifted with t\cm 
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FIG. 4: (Color online) Pseudorapidity distribution of charged 
hadrons in the FTPC acceptance for 0-20%, 20-40%, 40-100% 
central d+Au events scaled with Npart (open triangles), p+p 
measurements [5^] unshifted (triangles) and shifted by rjcM 
(lines) are overlaid (for further details see text). 



Centrality class 


tjcm 


0-20% 


-1.29 


20-40% 


-0.98 


40-60% 


-0.68 


60-80% 


-0.37 


80-100% 


-0.14 


40-100% 


-0.48 


0-100% 


-0.92 



TABLE II: t\cm defined as the weighted mean of the PHO- 
BOS dN ch /drj distribution for various centrality classes [24T [. 



tion: 



VCM 



) and the unshifted dN^/dr) distribu- 



SdAu(r],T]CM) 



dN%/d V \ 



V—Vcm 



dN%/dn\, 



(V, Vcm) ■ (2) 



Our ansatz is then that the suppression and enhance- 
ment of RdAu (and Rc p ) is mainly caused by this geo- 
metric asymmetry. It would follow from this ansatz that 
the difference in the observed centrality dependence of 
RdAu should be accounted for by simply scaling with 
SdAuiViJJCM) (Eq. The BRAHMS R dA u measure- 
ments [1J| at different pseudorapidities (t?=0, 1, 2.2 and 
3.2) on the d-side are then consistent with a universal be- 
haviour, reaching binary scaling at p t > 3 GeV without 
a significant Cronin enhancement at intermediate pt- 

A similar procedure can be applied to describe Rcp 
measurements for different centrality classes and pseudo- 
rapidities on the d- and Au-side of a d+Au collision. One 
has to modify SdAuiViVCM) (Eq. 2) to take the asym- 
metry in particle production - still visible in peripheral 
d+Au collisions (see Fig. - into account. This is re- 
alized by taking the peripheral d+Au dN c h/drj distribu- 
tion as reference and shifting the denominator in Eq. [5] 
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FIG. 5: (Color online) Rcp of charged hadrons at 77 ~ —3.1 
(Au-side) and at 77 « 3.1 (d-side) as function of pt for different 
centrality classes scaled with Scp{ri,riCM) (Eq.[3]). Midrapid- 
ity scaled charged hadron Rcp (0-20%) (\r]\ < 0.5) is shown 
as stars [2l]]. Also plotted are scaled BRAHMS Rcp mea- 
surements (0-20%/60-80%) HH in different 77 intervalls. 



according to the peripheral rjcM value: 



Scp(r],VCM) 



dN p c p h /d V \ 
dN™/d V \ 



V-VCM 



■(ViVcm)- (3) 



Rcp measurements from the FTPCs at \rj\ m 3.1 for 
different centrality classes scaled with Scp(j],t]cm) are 
shown in Fig. [5j These results suggest that the observed 
enhancement of Rcp on the Au-side for different cen- 
trality classes can be explained by an enhancement of 
particle production caused by the nuclear stopping of 
the deuteron and described by a shift of the center-of- 
mass with respect to the most peripheral d+Au collisions. 
The suppression of Rcp on the d-side is explained anal- 
ogously. Universal scaling behaviour of Rcp scaled with 
ScpIv^Vcm) at \r)\ ~ 3.1 is visible for p t > 1 GeV/c in 
the FTPC measurements (see Fig. [5}. In addition scaled 
BRAHMS Rcp measurements for the top 20% central 
d+Au events at different pseudorapiditics are also over- 
layed in Fig. [5] Universal scaling behaviour of scaled 
Rcp is seen in both STAR and BRAHMS measurements 
reaching binary scaling for p t > 2 GeV/c. The deviation 
from the supposed scaling for the BRAHMS Rcp mea- 
surements at midrapidity (and also the discrepancy to the 
STAR midrapidty measurements) might be explained by 
the different centrality definitions, where a bias to higher 
particle multiplicities might be present in the BRAHMS 
data as discussed in section IIII CI It should be noted 
that systematic errors in the scaling due to uncertain- 
ties in the determination of t]cm were not estimated in 
this analysis. One would expect an overall normaliza- 
tion uncertainty, but the main feature of scaling should 
be preserved. The simple data driven picture outlined in 
this paper was meant to point out the importance of the 
collision geometry in d+Au collisions. 



V. CONCLUSIONS 

Centrality dependence of the dN/dr) distributions in 
v / saT/v=200 GeV d+Au collisions is presented. Within 
errors, dN/drj cannot discriminate between different 
model calculations. However, our data show an increase 
of (pt) as a function of centrality on the Au-side, whereas 
on the d-side no centrality dependence is visible. This 
result together with the suppression of Rcp on the d- 
side and its enhancement on the Au-side relative to 
mid-rapidity cannot be described consistently by current 
model calculations [TBI ]. A similar study of comparing 
particle production at high p t in d- and Au-side at mid- 
rapidity [29l | has ruled out models based on incoherent 
initial multiple partonic scattering and independent frag- 
mentation. It also showed that models based on nuclear 
shadowing incorporating extremes of gluon shadowing at 
low x can not account for the difference in particle pro- 
duction in d- and Au-side at high p t at mid-rapidity. 

This paper demonstrated, that in a simple stopping 
picture the main features of the pscudorapidity and cen- 
trality dependence of Rcp (and RdAu) at higher rapidi- 
ties can be explained by the suppression (enhancement) 
of particle yields in d+Au relative to p+p collisions and 
peripheral d+Au collisions. Simulation studies show a 
small effect of shadowing in HIJING on Rcp, especially 
on the d-side, supporting this geometric picture. This re- 
sult is also confirmed by measurements of R p pb by NA49 
[30l | at an order of magnitude lower energies, where shad- 
owing and gluon saturation are expected to be small, 
which show qualitatively the same characteristic central- 
ity and rapidity (xp) dependence and are consistent with 
the stopping picture. 

Taking the asymmetry in particle production 
characterized by a shift of the center-of-mass in the 
asymmetric d+Au collisions with respect to the nucleon- 
nucleon or peripheral d+Au center-of-mass system - 
into account, Rcp (and RdAu) show a universal scaling 
behaviour independent of centrality and pseudorapidity 
in d+Au collisions at RHIC energies. On the other 
hand, the success of the CGC saturation approach 
[3l| and stopping picture in quantitatively describing 
BRAHMS Rcp (and RdAu) measurements could be 
interpreted as a link between saturation as the possible 
origin of nuclear stopping. In that case deviations from 
the observed empirical scaling in asymmetric collision 
systems at different energies can help to quantify the 
onset of saturation effects in heavy-ion collisions. 
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